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Abstract. We expose a relationship between jamming and a generalization 
of Tutte's barycentric embedding. This provides a basis for the systematic 
treatment of jamming and maximal packing problems on two-dimensional 
surfaces. 



1. Introduction 

In a seminal paper Q], W. T. Tutte addressed the problem of how to embed a 
three-connected planar graph in the plane. He proposed to fix the positions 
of the vertices on one face and to let the other (inner) vertices be barycenters 
of their neighbors (see Fig. 0. The barycentric embedding is unique. It 




Fig. 1. Barycentric embedding of a graph with TV = 13 vertices (three outer vertices). 
minimizes the energy 



* = Ei 



i<j 



over the positions r j of the inner vertices. 

The purpose of this paper is to study jamming, which is of importance 
for the physics of granular materials and of glasses EII5E1 . and has many 
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applications in mathematics and computer science 0. We expose a re- 
lationship between jamming and a generalization of the barycentric em- 
bedding, and provide a basis for the systematic treatment of jamming on 
two-dimensional surfaces. 

Ensembles of non-overlapping disks of equal radius may contain sub- 
ensembles of disks which do not allow any small moves, regardless of the 
positions of the other disks: In Fig. |2] (showing disks in a square), disks 
i, j, k, I, m, n are jammed, while o is free to move. 




3t 

>i A— it 

ii ip — 



Fig. 2. Left: Configuration of seven disks in a square. Disks i,j,k,l,m,n are jammed. 
Right: Contact graph of the jammed sub-ensemble. Edges among outer vertices are omitted. 



The position r-j of the center of disk i must be more than a disk diam- 
eter away from other disk centers, and more than a disk radius from the 
boundary. In a jammed configuration, r-j locally maximizes the minimum 
distances to all other disks, and twice the distances to the boundaries. For a 
disk i not in contact with the boundary, we have 



mm r,- 



max mm r — i"j 



Equivalently, r j is a local minimum of the repulsive energy 



£(r) 






..19 



(1) 



in the limit q — » oo. Empirical approaches with a repulsive energy, as in 
eq. (Q, with very large q have been important to actually find jammed con- 
figurations [6|. As the minimum is local, one cannot prove that with this 
method all jammed configurations are generated. 

The relationship between jamming and the geometric representation of 
graphs was first pointed out, a long time ago, by Schiitte and van der Waer- 
den [7|. Each jammed disk corresponds to an inner vertex of a graph, and 
each touching point with the boundary to an outer vertex. The edges of the 
graph refer to contact of disks among themselves and with the boundary, as 
shown in Fig. 12 

In this paper, we show that the position r j is not only the local maximum 
of the minimum distance to all other disks, but also the global minimum of 
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the maximum distance to all the neighbors of i 

max |r.j— r a |=min max |r — r a |. 

a=j,k,l,m r a=j,k,l,m 

This leads us to study M. -representations of graphs (as in Fig. with 
inner and possibly outer vertices, where each one minimizes the maximum 
(rescaled) neighbor distance rather than the mean squared distance, as in 
Tutte's barycentric embedding. Outer vertices are either fixed or restricted 
to line segments. 




Fig. 3. Stable M -representation of the graph of Fig.Q with identical positions of the outer 
vertices. Three faces of this representation are flat. 



We define stable representations as .M-representations which are lo- 
cal minima with respect to an ordering relation. This relation replaces the 
notion of an energy which cannot be defined in this setting. 

We establish that the problem of finding a stable representation in the 
plane, torus, or on the hemisphere has an essentially unique solution for 
any graph. We show that M. -representations of three-connected planar and 
toroidal graphs are convex pseudo-embeddings, and that the set of regular 
three-connected stable representations contains all jammed configurations. 
This puts jamming in direct analogy with the barycentric embeddings. On 
the sphere, stable representations are not unique, but we conjecture that 
their structure is restricted. 

One application of jamming is the generation of packings of N non- 
overlapping disks with maximum radius. Such a maximal packing con- 
tains a non- trivial jammed sub-ensemble, since otherwise we could increase 
the radius of each disk. The remaining disks of a maximal packing are 
not jammed (as disk o in Fig. |5J, and confined to holes in the jammed 
sub-ensemble. In these holes, we can again search for jammed configura- 
tions with suitably rescaled radii. This gives a recursive procedure to com- 
pute maximal disk packings, which relies on the enumeration of (three- 
connected) planar or toroidal graphs and a computation of their jammed 
representations which form a subset of the stable representations. Practi- 
cally, we generate the stable representations with a variant of the minover 
algorithm [ 8 1 which appears to always converge to a stable solution, on the 
plane, torus, and on the sphere. 

The most notorious instance of maximal packing is the N = 13 spheres 
problem for disks on the sphere. It has been known since the work of 
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Schiitte and van der Waerden [7] that 13 unit spheres cannot be packed 
onto the surface of the unit central sphere (a popular description has ap- 
peared recently in the French edition of Scientific American [9 1). However, 
the minimum radius of the central sphere admitting such a packing is still 
unknown, as it is for all larger N, with the exception of N = 24 [10]. The 
problem of packing spheres on a central sphere is clearly equivalent to the 
problem of packing disks on a sphere. 

Our strategy for solving the maximal packing problem will be complete 
once the following conjectures are validated: 

Conjecture 1. There exists a finite algorithm to find a stable representation 
of a given graph. 

Conjecture 2. Each graph on the sphere with a fixed set of edges crossing a 
given equator has at most one non-trivial stable representation up to sym- 
metry transformations on the sphere. Furthermore, jammed configurations 
are stable. 

At present, we are able to prove Conjecture |2] for a fixed representation 
of edges across an equator, rather than their set. The conjecture is backed 
by extensive computational experiments. For planar region and torus, only 
Conjecture [I] is needed. 



2. .M-representations 

In this section we discuss representations of graphs in a planar region, on 
the torus, and the sphere. By torus we mean a rectangular planar region 
where the parallel sides are formally identified. In a representation, each 
vertex is a point, and each edge the shortest connection between vertices. 
On the torus, the rectangle can always be chosen so that vertices do not 
lie on its sides. We require that in a representation the shortest connection 
between vertices connected by an edge is uniquely determined. A represen- 
tation is an embedding if it corresponds to a proper drawing. We recall that 
a graph is fe-connected if it has more than k vertices and remains connected 
after deletion of any subset of k — 1 vertices. Furthermore, we use two basic 
facts of graph-theory: the faces of a two-connected planar embedding are 
bounded by cycles, and embeddings of a three-connected planar graph have 
a unique list of faces and incidence relations. 

Each toroidal representation of a graph gives rise to a unique periodic 
representation by tiling the plane with the rectangles, as shown in Fig.HJ A 
proper toroidal representation has edges crossing each side of the rectangle 
and no outer vertices. 

Definition 1 (Inner and outer vertices). Each graph contains a possibly 
empty subset O of outer vertices. Each outer vertex is fixed or constrained 
to lie on a line segment such that any choice of outer vertices forms a sub- 
division of a convex n-gon, n < \0\. 
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Fig. 4. The periodic representation of a toroidal graph with six vertices. 

As indicated in the introduction, we define a rescaled distance, in order 
to treat outer and inner vertices on the same footing. 

Definition 2 (Rescaled distance). The distance between vertices i and j is 



L JI7 



lij 



r, 



L Jh 



where 7^ are positive constants, and | | denotes the Euclidean distance. 
In a given representation, we denote by 1(e) the (rescaled) length of an edge 
e, i.e. the distance between its end-vertices. 




Fig. 5. Representation of planar graph in the plane. The outer vertex i is constrained to lie 
on a line segment, whereas j and k are fixed. 



2.1. A4- center of vectors 

Definition 3 (Radius). The M— center of a finite number of vectors r^, i € 
/ is the vector r^ minimizing the radius ofr: p(r) = maxj |r — rj| 7 



p(roo) = min p(r). 



Note that the .M-center is uniquely determined: If there were two M.- 
centers with the same radius p, then the intersection of the corresponding 
circles of radius p would contain all the neighbors, but this intersection is 
contained in a circle of smaller radius. 
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Lemma 1 (No local minimum besides global one). If r is not the M.- 

center of vectors n, i £ I, then for each 5 > there is a vector r' with 
|r' — r| < 5 such that p(r) > p(r'). 

To determine the .M-center of vectors r , , i G /, we may consider all pairs 
and triples of vectors, construct the (unique) circle passing through them 
and check that no vertex is outside the circle. The center of the smallest 
circle such obtained is the M. -center. 

Definition 4 (.M-representation). An Ai-representation of a graph is a 
representation where each inner vertex is the M— center of its neighbors. 

Definition 5 (Pseudo-embedding). A representation of a graph is called 
pseudo-embedding if it is an embedding except that some faces may col- 
lapse into a line. Such faces will be called flat Moreover, a convex pseudo- 
embedding has convex faces and each flat face is a topological subdivision 
ofCi, a cycle of length two. 

An example of a convex pseudo-embedding is shown in Fig. |3] 

Proposition 1. Let £ be a representation of a three-connected planar graph 
on a plane or hemisphere, such that each inner vertex belongs to the con- 
vex hull of its neighbors, with non-empty set of outer vertices. Then £ is a 
convex pseudo-embedding. 

Proof. We proceed analogously to the paragraphs 6-9 of [1]. Since G is 
three-connected, its set of faces is uniquely determined, and each face is 
bounded by a cycle. O denotes the set of outer vertices. 

Let I be a line in the plane or a non-trivial intersection of a plane with 
the hemisphere and define g(v),v £ V, as the perpendicular distance of v 
to I, counted positive on one side and negative on the other side of I. 

The outer vertices with the greatest value of g are called positive poles 
and those with the least value of g are negative poles. Note that the sets of 
positive and negative poles are disjoint since O forms a subdivision of a 
convex n-gon. 

A simple path P = vi, ■ ■ ■ ,-Ufc of G is right (left) rising if for each i, 
g{vi) < g(vi+i) or g(vi) = g(v i+1 ) and v i+ i is on the right (left) hand-side 
of V{. Right (left) falling paths are defined analogously. 

Lemma 2. Each vertex v of G different from a pole has two neighbors v' 
and v" so that g(v') < g(v) < g[y") or g(v') = g(v) = g[y") and v 
belongs to the line between v' and v". 

Proof. This follows for outer vertices since they form a subdivision of a 
convex n-gon, and for inner vertices because of the convexity assumption. 

Lemma 3. Let v be a vertex of G. There is a right rising and a left rising 
path from v to a positive pole, and also both right and left falling paths 
from v to a negative pole. 
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Proof. By LemmalJ]?; has a neighbor v' with g(v') > g(v) or g(v') = g{v) 
and v' is on the right hand-side of v. Since G is three-connected, v' has a 
neighbor different from v. Using Lemma|2l we can monotonically continue 
from v' . This constructs a right rising path, and the remaining paths may be 
obtained analogously. 

Lemma 4.1fv ^ O then v belongs to the convex hull ofO. 

Proof. If such v does not belong to the convex hull of O, then let / be a line 
in the plane (cycle on the hemisphere) which defines a separating plane, 
and we get a contradiction with Lemma |3] 

Lemma 5. Let F be a face of G and vi,v[,V2,v' 2 vertices of F appearing 
along F in this order. Then G does not have two disjoint v% , v 2 and v[ , v' 2 
paths. 

Proof. This is a simple property of a face of a planar graph. 

Lemma 6. If a face F is flat then it is a topological subdivision of C 2 . 
Furthermore, let e be an edge of a face F and let I be a line in the plane (a 
cycle on the hemisphere) containing e. Then F is embedded on one side of 
I. 

Proof. This simply follows from Lemma|3]and Lemma|5](see Fig. |6j. 



j^d 




Fig. 6. Left: a flat face must be a subdivision of Ci, Right: each face must lie on one side 
of incident edge e. 



It follows from Lemma that each face is a subdivision of a convex 
n-gon or flat and subdivision of Ci- 

Each edge belongs to exactly two different faces. An edge is redundant 
if it belongs to two flat faces. More generally, in a two-connected represen- 
tation with prescribed faces such that each flat face is a subdivision of Ci, 
a path which is a subdivision of an edge is redundant if it belongs to two 
flat faces. A graph is a simplification of G if some redundant edges and, 
thereafter, maximal redundant paths have been deleted. 

Lemma 7. Aflat face of a simplification of G is a subdivision ofC 2 . 

Proof. If we delete e and unify the two faces containing e, we get a pla- 
nar graph. If the statement does not hold then we can again use Lemma |3] 
and Lemma |5] to obtain a contradiction. The same applies for a maximal 
redundant path. 
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Let G' be the smallest simplification of G and let F be a flat face of G'. 
We know that it is a subdivision of C2, and each edge of F belongs to one 
of the two sides of C2. 

Lemma 8. Let e be an edge of G' and let I be the line in the plane (cycle 
on the hemisphere) containing e. 

1. Ife belongs to aflat face F, then the faces incident with edges of differ- 
ent sides of F are on opposite sides of I. 

2. If edge e does not belong to aflat face, then the two faces incident with 
e lie on opposite sides of I. 

Proof. For the second property: as in the proof of Lemma if we delete 
e and 'unify' the two faces containing e, we get a planar graph. If the two 
faces lie on the same side of e, we can use Lemma |5] The first property is 
analogous. 

Let |Gr| denote the subset of the surface consisting of the embeddings of 
the vertices and edges of G, and let S denote the complement of |G|. 

We define a function d on S as follows: d(x) = 1 if x is not within 
the convex hull of O, otherwise, d(x) equals the number of interiors of 
faces to which x belongs. The correctness of this definition is guaranteed 
by Lemma|4] 

Lemma 9. For each x G S, d(x) = 1. 

Proof. It follows from Lemma[8]that the function d does not change when 
passing an edge. However, it cannot change elsewhere and outside of the 
convex hull of O it equals to 1. Hence it is 1 everywhere. 

Lemma 10. If an edge e intersects the interior of an edge e', then one of 
them is not in G' or they belong to opposite sides of aflat face of G 1 . 

Proof. This is a corollary of Lemma|9| 

□ 

Corollary 1 (.M-rep. is pseudo-embedding). An M.-representation with- 
out outer vertices of three-connected planar graphs on a sphere or three- 
connected proper toroidal graphs on a torus is a convex pseudo-embedding. 

Proof. As the number of vertices is finite, we can always find a cut (rect- 
angle and plane through the center, respectively), which does not contain 
any intersection of two edges. The corollary follows by taking as outer ver- 
tices the intersection of edges with the cut. If the cut does not intersect any 
edges, the representation is trivial. 

Definition 6 (Ordering of representations). Consider two representations 
£ and £' of a graph G We say that £ is smaller than £' (£ < £') if the 
ordered vector of lengths of the edges containing an inner vertex of £ is 
lexicographically smaller than the ordered vector of lengths of the same 
edges in £'. 
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The above ordering relation cannot generally be mapped into the real 
numbers, because the real axis does not admit an uncountable number of 
disjoint intervals. Therefore, there is no 'energy' (generalizing eq. Q) such 
that£ <£' ^E(£) < E{£'). 

Definition 7 (Stable representation). Consider a representation £ of a 
graph G = (V, E) with inner, and possibly outer vertices i at positions 
Ti. £ is stable if there exists a value 5 such that all embeddings £' ofG with 
vertices at r£ with |rj — r^| < 5\/i satisfy £' > £. 

Proposition 2. Stable representations are M— representations. 

Proof. Let the vertex i of £ have the radius pi and let edge {i, j} have 
length pi. Note that pj > pi. If i is not the .M-center of its neighbors, then 
it follows from Lemma [2 that there is a representation £' obtained from £ 
by a small move of vertex i, such that p\ < p^. All edges {k, 1} with length 
bigger than pi are the same in £ and £'. No edge {k, 1} of length pi in £ is 
longer in £' and at least one such edge has shortened. Finally, edges {k, 1} 
shorter than pi in £ may become longer in £'. As a result, we have £' < £, 
which is impossible for a stable representation. □ 

Proposition 3 (Existence of stable representation). Each graph has a sta- 
ble representation. 

Proof. Let 5 > be a sufficiently small constant. Define a sequence of rep- 
resentations <?i, <?2 ; • • • as follows: £\ is arbitrary. If £i is unstable let £j + i 
be a lexicographically minimal representation where each vertex has moved 
by at most 5 (it exists by compactness). In particular £,- i+ i < £i. Again by 
compactness, there is a converging subsequence of representations £'■ with 
limit £'. £' must be stable since otherwise for £[ very near to £', there is 
a closeby representation £ < £'. Taking into account the minimality rule 
in the construction of the sequence of representations, this contradicts the 
assumption that £[ monotonically decreases in lexicographic order to £' . □ 



2.2. Uniqueness of stable representations 

Proposition ffl implies that each .M-representation is a convex pseudo- 
embedding. Whereas Tutte's barycentric embedding is unique, the M.- 
representations are not necessarily unique, as can be seen by the counter- 
example sketched in Fig. |5] (the vertices of the inner triangle are in Ap- 
position; they can be rotated and rescaled, to remain in M. -position). How- 
ever, there is a unique stable representation. 

Lemma 11. Consider two-dimensional vectors ri , r2, r' l5 r' 2 with 

ri = (xi,yi), etc 
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and two midpoints f i and r 2 : 

xi = \{xi +x[) 

Vi = \(yi + v'i)- 

We then have 

|ri -r 2 | 7 + |r' 1 -r 2 | 7 



|ri -r 2 | 7 < 



2 



We have \t\ — r 2 | 7 = |r^ — r 2 | 7 = |ri — r 2 | 7 only for parallel transport: 
t\ = r[ + c; r 2 = r 2 + c. 

Proof. Follows from triangle inequality |a+b| < |a| + |b|, witha = ri— r 2 
and b = r[ — r 2 with equality only for parallel transport. 

Note that if \r% — r 2 | / |r^ — r 2 |, the midpoint distance |ri — r 2 | is smaller 



Proposition 4 (Unique stable representation in the plane). Each graph 
G has a unique stable representation in the plane (up to parallel transport). 

Proof. We assume the contrary. Let representations £° and £ l , realized by 
vectors r? and r\, be two stable representation. We can assume £° < £ l . 
Consider the representations £ a realized by 



r 



+ ax[r J 1 -r°] < a < 1. 



The representations £ a exist. We denote by e° and e 1 the representations 
of edge e in £° and £ , respectively. Let e\, . . . , e^ be the ordered vector 
of edge lengths. Let k be the smallest index such that e\ is not parallely 
transported to e\. We observe the following: if e is an edge of G such that 
/(e 1 ) = l( e \), then Z(e°) < l{e\) since £° < £ l . It means by Lemma [TT1 
that £ a < £ x Va < 1, which implies that £ x is not stable. □ 

Proposition 5 (Unique stable representation on torus). Each graph G 
has a unique stable representation on the torus if the sets of edges crossing 
each boundary are prescribed {up to parallel transport). 

Proof. The representations £ a of the previous proof can analogously be 
applied to the corresponding periodic representations, both for edges in the 
inside of one rectangle and for the edges going across the boundary. □ 

This means that the number of stable representations of a toroidal graph is 
bounded by the number of possible sets of boundary horizontal and vertical 
edges. 

Next we will discuss uniqueness of stable embeddings on the hemi- 
sphere. The following Lemma[T2lis a nontrivial variant of Lemma^J Ob- 
viously, the triangle inequality remains valid in three dimensions, but the 
midpoint would not lie on the surface of the sphere. 
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Lemma 12. Consider three-dimensional vectors ri,r2,ri,r 2 on the unit 
hemisphere with 



ri = (>i,yi,zi = +yl -x{ -y{), etc. 

and two midpoints r i and Y 2 which are defined with respect to a projection 
of vectors on the equator z = 

Xi = \{xi + z-); y i = \{y { + y-); Zj = yjl - xf - yj i = 1, 2. (2) 



We ?/2£?i /zave, for the three-dimensional Euclidean distance-squared: 

1- - ,2 ^ |ri - r 2 | 2 + [r; - r 2 | 2 
r i- r 2 < ?; • 



(3) 



In @, we can have equality only for generalized parallel transport with 
x\ — x[ = c(x2 — x' 2 ) and y\ — y[ = c(y2 — y^far special values of c. 



Proof. We can write © as 



{-Zl-^2? < 



x\ + x[ x 2 + x' 2 



(gi - x 2 ) 2 Qi - x' 2 ) 2 



+ same terms in y, y' + \{z\ — Z2) 2 + \(z[ — z' 2 ) 2 . (4) 
Explicit calculation shows that the terms on the first row of expression © 



is 



Xi + x'i X 2 + x' 2 



, {x 1 -x 2 )' 2 | K-x 2 ) 2 



= \[X\ -X 2 ~X[ +X 2 } 2 , 

which allows to show that © and (J3} w& equivalent to 

(zi - z 2 ) 2 <\{xi-x 2 - x[ + x' 2 ) + \ (yi - 2/2 - y'i + 2/2) 

+ \{z l -z 2 ) 2 + \{z' 1 -z' 2 )\ (5) 

Furthermore, we have, from eq. (|2|i 



z 2 = l 



»i + x'i \ fyi + yl 



2 / V 2 
Uxt - x'^ 2 + \{ Vi - y'^ 2 + i(z 2 + z' 2 ) i = 1, 2. 
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The inequality © now follows from the triangle inequality (|a| — |b|) 2 < 
| a — b| 2 with the four-dimensional vectors 

z\ z[ x\ - x[ yx -y[ 



s/2 ^ 2 ' 2 

/^2_ _4_ ^2 ~ 4 2/2 -2/2 

VV2'V2' 2 ' 2 

which evidently satisfy |a| = zx, |b| = Z2 and |a — b| 2 equal to the r.h.s. 
of ©. □ 

Proposition 6 (Unique stable rep. on hemisphere). Each graph G on the 
hemisphere with fixed outer vertices has a unique stable representation. 

Proof. Using the definition of midpoints on the sphere from Lemma[l2l we 
can define valid representations £ a as in the proof of Proposition HJ Then 
we still have £ a < £ l \/a < 1. It is easy to see that, with fixed outer 
vertices, generalized parallel transport is impossible. □ 



2.3. Stable representations on the sphere 

On the sphere, there can be several non-trivial stable representations. To see 
this, consider the equator-representation of Fig.0 Besides a central cycle, 



Vj 




Fig. 7. An equator embedding of a graph G, with a central cycle at z — 0, and upper and 
lower subgraphs G upp er (at z > 0) and Glower (at z < 0). The edges on the central cycle are 
longer than all other edges. 



at z = 0, there are vertices in the upper subgraph G U pper (with z > 0) and 
in the lower subgraph Gi OW er (at z < 0). Furthermore, we suppose that the 
edges on the central cycle are longer than those in the rest of the graph. 
An equator representation can give rise to two inequivalent representations, 
namely by pulling the central cycle up to z > 0, or down to z < 0. 

Proposition |6] allows us to observe that nevertheless, some degree of 
uniqueness can be preserved. 

Proposition 7 (Unique representation with fixed cut). There is unique 
stable representation of a graph on a sphere when the edges crossing a 
given equator are fixed. 
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As mentioned in the introduction, extensive computing experiments sug- 
gest the stronger statement of Conjecture 13 This would imply that a given 
graph has only a finite number of stable representations (up to symmetry 
operations). 



3. Jamming 

In this section we derive basic properties of jammed configurations of disks 
on planar region, torus and sphere. 

Definition 8 (Jammed embedding). An embedding £ of a graph is jammed 
if 

1. £ belongs to the convex hull of the set O of outer vertices, and each 
outer vertex is connected to at least one inner vertex. 

2. £ is regular, i.e. all edges have length 7^ x d and the distance between 
any two vertices k, I not connected by an edge is strictly bigger than 
7kix d. 

3. there is a 5 such that no representation £' with | r^ — r , | < 5 has some 
edge longer and no edge shorter than in £. 

The three central disks in Fig. [8] are not jammed, even though each one 
cannot move individually. 




Fig. 8. Left: Configuration of five disks, in which no disk can move by itself, but three disks 
can move together, as indicated. Right: The (unjammed yet stable) embedding correspond- 
ing to the configuration. Edges among outer vertices are omitted. 



Let us recall that in a representation, set O of outer vertices forms a 
subdivision of a convex n-gon. Let us denote this cycle of outer vertices by 
Cq and if G is a jammed embedding then let us denote by Go the pseudo- 
embedding obtained from G by adding the edges of Co', let us note that 
Co bounds the outer face of Go- In a jammed embedding G each inner 
vertex has degree bigger than two and at most five on the sphere, and at 
most six on planar region and torus. If Go is two-connected then each 
inner face is convex ([7]). Below we show that each jammed graph Go is 
three-connected. 

Lemma 13. If a graph G is connected but not two-connected and has no 
vertex of degree 1, then there are two vertices v\,v 2 (possibly v\ = v 2 ) and 
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components GiofG\vi(i = 1, 2) such that each d U Vi is two-connected 
and G\ U v\ is a subgraph ofG\ G2. 

Lemma 14. If a graph G is two-connected but not three-connected, and 
has no vertex of degree 2, then it has two pairs Vi, V2 of vertices (possibly 
with Vi fl V% 7^ 0) and components G{ofG\Vi(i = 1, 2) such that G L U V- L 
is two-connected and G\ U V\ is a subgraph ofG\ GV 

Definition 9. Let graph G be a two-connected regularly embedded graph, 
let F be a face of G bounded by a cycle C and let v £ C. We say that v is 
convex with respect to F if there is a plane containing v and perpendicular 
to the surface (planar region, torus, sphere), so that a small neighbourhood 
ofv in F lies completely in one half -space defined by the plane. 

We note that if F is a convex region bounded by a cycle then each vertex 
of the cycle is convex with respect to F. 

Lemma 15. Let F\ , F2 be connected regions bounded by cycles C\ , C2 and 
such that none of them covers the whole planar region (torus, sphere, hemi- 
sphere), but F\ U F2 do. Then there are at least three non-convex vertices 
of C\ with respect to F\ or of C<i with respect to Fi, 

Proof. Cycle C2 must be embedded inside F\ and F2 must be the region 
defined by C2 that is not a subset of F\ . Then C2 is a non-self-intersecting 
cycle on F\ . As such, it must have at least three sharp corners on F\ . 

Proposition 8 (Jammed graph three-connected). If an embedding G is 
jammed, then Gq is three-connected. 

Proof. Let Go be a minimum counter-example. If Go is not connected 
then each component is jammed and three-connected by the minimality 
assumption. Consider the embeddings of different components G\ and G2 
in the embedding of G. G\ is completely embedded in one of the faces 
of G2 and vice versa, which is not possible by convexity of faces and by 
Lemma[T5l 

If Go has a vertex of degree at most two then it cannot be jammed. 
Therefore, we suppose that Go is without a vertex of degree two and either 
connected or two-connected. By Lemma[T3land LemmaElthere is a subset 
of vertices V\ and a component G\ of Go \ Vi such that G\ U V\ is two- 
connected and has one of the following two properties: 
l.Vi consists of a single vertex v\ and there is a vertex V2 with V2 = {^2} 
and a component G2 of Go \ V2 so that G2UV2 is a two-connected subgraph 
of Go\Gi. 

2. V\ consists of two vertices and there is a subset V2 of two vertices and a 
component G2 of Go \ V2 so that G2 U V2 is a two-connected subgraph of 
Go \ G\. 

We consider the embedding of G\ U V\ induced by the embedding of 
Go (see Fig.©. Let F\ be the face in which Go \ G\ is embedded and let 
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C\ be the bounding cycle of F\ . Clearly V\ C C\ . Moeover F\ cannot be 
the outer face of the embedding of Go since Go \ G\ is embedded there. 
Hence all the vertices of C\ \ V\ are convex with respect to F\ . 

Let F2 be the face of the induced embedding of G 2 U V 2 which contains 
C\ and let C2 be the cycle that bounds F 2 . Clearly V 2 C C 2 and as above, 
F 2 cannot be the outer face of the embedding of Go- Hence all vertices of 
C 2 \ V2 are convex with respect to F 2 . Then F\ , F% , C\ , C 2 satisfy the prop- 
erties of Lemma [131 but V\ and V2 have only two vertices, a contradiction. 
D 




¥ 



:v, 



]^2 




Fig. 9. 



Proposition 9 (Jammed graphs stable). Jammed embeddings are stable 
on the planar region and torus and on a hemisphere. On the sphere, jammed 
embeddings are stable if we fix representations of edges across an equator. 

Proof. We show that an unstable regular embedding £ has a small move 
which leaves some edges the same and increases the others, which means 
that it is not jammed. For any 8 > 0, there is a representation £' within 5 
of £ such that some edges decrease in length, and the others stay the same. 
Let £" be the inverse of the move which took £ into £' (on the hemisphere: 
inverse of the projected move). The inverse move exists, since a jammed 
embedding in the plane or on the torus cannot have an outer vertex on an 
end point of the corresponding line segment, and for the hemisphere, if 
originally an inner vertex positioned on the equator moved, then £ could 
not have been jammed. For the representation e of an edge in £, let e' and 
e" be the respective representations of the same edge in £' and £". We 
can use Lemma ^2 (on planar region and torus) and Lemma \^\ (for the 
hemisphere) to show that l(e') = 1(e) =^> l(e") > l( e ) an d K e> ) < 
l( e ) =^ l( e ") > /( e ). □ 

Note that the converse of Proposition |9]is not true and that stable represen- 
tations are not necessarily jammed. An example is shown in Fig. [8] 

As mentioned in the introduction, extensive computing experiments sug- 
gest, for the sphere, the stronger statement of Conjecture [21 
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4. Algorithms 

In Section ITT1 we discussed a finite algorithm for the determination of the 
.M-center of vectors r-j. This algorithm is of practical use because a circle 
in d-dimensional space is already specified by d+1 points. The incremental 
minover algorithm [ 8 1 remains useful in high dimension d and is trivial to 
implement. For a finite number of vectors r j , i G / on the unit sphere, it is 
defined by 

Ro = 0, Rfc+i <— R/c + r-j m . n , 

where the index i m j n is a minimal overlap (scalar product) vector with 



(Rt.r. . 



min(R fc ,rj 



It can be proven 1 8 ] that 

Rfc/|Rfc| ^roo 

under the condition that a vector R exists with (R, r-j) > \/i G /. 





Fig. 10. Left: Minimum overlap algorithm. The rescaled vector Rfc/|Rfc | converges to the 
.M-center of vectors on the sphere. Right: On the plane, the move is in direction of the 
vector r, with maximum distance to Rfe. The amplitude e^ of the move decreases with k, 
but J2k £ fc diverges. 



In the minover algorithm (on the sphere), the corrections to R^ decrease 
with increasing k. On the plane or the torus, we can do the same by using 
an update 

Rfe+i <— Rfc + £fc [r« max - Rfc], 

where rj max is the vector of maximum distance to R^ (see Fig. [lOl . and 
where the sequence e^ satisfies the conditions: 

efc —>■ for k —>■ cxd 

y efc — > oo for k — y cxd. 
k 

This algorithm was applied to all vertices sequentially in order to compute 
stable M. -representations. 
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As a simple test, we have run this algorithm on the three-connected 
graph of Fig. [2 and Fig. |3] embedded on a sphere and starting from an 
equator position. This graph, incidentally, corresponds to the conjectured 
optimal packing for the thirteen-sphere problem. The algorithm converges 
rapidly to the conjectured optimum solution [6], which is shown in Fig. [HI 




Fig. 11. Left: Conjectured optimal configuration of 13 disks on a unit sphere, and corre- 
sponding representation of vertices, obtained by computational experiment as stable M- 
representation of the graph of Fig.rjjand Fig.|5] The algorithm of Section|4]was used. 



As stated in Conjecture^ we are convinced that a finite algorithm for 
computing a stable J\A -representation exists. 
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